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ABSTRACT 

Among the initial results from Kepler were two striking light curves, for KOI 74 and KOI 81, in which 
the relative depths of the primary and secondary eclipses showed that the more compact, less luminous object 
was hotter than its stellar host. That result became particularly intriguing because a substellar mass had been 
derived for the secondary in KOI 74, which would make the high temperature challenging to explain; in KOI 8 1 , 
the mass range for the companion was also reported to be consistent with a substellar object. We re-analyze 
the Kepler data and demonstrate that both companions are likely to be white dwarfs. We also find that the 
photometric data for KOI 74 show a modulation in brightness as the more luminous star orbits, due to Doppler 
boosting. The magnitude of the effect is sufficiently large that we can use it to infer a radial velocity amplitude 
accurate to Ikms"'. As far as we are aware, this is the first time a radial- velocity curve has been measured 
photometrically. Combining our velocity amplitude with the inclination and primary mass derived from the 
eclipses and primary spectral type, we infer a secondary mass of 0.22 ±0.03Mq. We use our estimates to 
consider the likely evolutionary paths and mass-transfer episodes of these binary systems. 

Subject headings: binaries: eclipsing — stars: evolution — stars: individual (KOI 74, KOI 81) — techniques: 
photometric — white dwarfs 



1. INTRODUCTION 

It has long been realized that precision photometry from 
space could allow detailed probes of stellar interiors through 
pulsations, and of planetary and stellar companions through 
transits, eclipses, and flux variations due to tides and irradi- 
ation. With the launch of Kepler, measurements with long- 
term precision of at least 1 part in 10'' are now being collected 
for an unprecede nted sample of stars (Caldwell et al. 2010; 
iKoch et al.ll2()T0h . 

Among Kepler'?, initial discoveries are two "objects of in- 
terest," KOI 74 and KOI 81, that have striking l ight curves, 
showin g both pr imary and secondary eclipses (iRowe et"an 
l2010al hereafter iRlOt where needed, we refer to the origi- 
nal version posted on arXiv as iRlOvll) . In both systems, the 
relative eclipse depths show that the more compact, less lu- 
minous object is hotter t han its stellar host. This result be- 
came puzzling as iRlOvll derived a clearly substellar mass for 
the secondary in KOI 74; their mass range for the compan- 
ion in KOI 81 was also consistent with a substellar object. If 
the smaller, hotter companions were indeed planets or brown 
dwarfs, then it would not be obvious how their properties 
could reasonably be explained. 

In this work, we re-analyze the Kepler data for KOI 74 and 
KOI 8 1 , and show that the photometric data for KOI 74 con- 
tain clear evidence for Doppler boosting with orbital phase, 
in addition to ellipsoidal light variations. We use this to make 
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what is, as far as we are aware, the first measurement of a 
radial-velocity amplitude from photometry. With our am- 
plitude, it is clear that the companion is not substellar, but 
rather has a mass of ^0.2Mo. For KOI 81, we find that the 
mass is likely similar. We demonstrate that all properties are 
consistent with the companions being low-mass white dwarfs 
(WDs), and discuss possible evolutionary histories for these 
systems. 

2. THE LIGHT CURVES REVISITED 

We retrieved the light curves of KOI 74 and KOI 8 1 from 
the archive, and show the raw aperture fluxes in Figure[T] For 
both sources, one sees the eclipses and transits superposed on 
roughly quadratic trends; in KOI 81, one also sees the pul- 
sations reported by lRlOl We first review what one can infer 
from the transits and eclipses, and then discuss additional in- 
formation seen at other phases. 

2.1. Constraints from the Transit and Eclipse Light Curves 

The transit an d ecl ipse light curves were fitted with de- 
tailed models by 'RlO, and we use their ephemerides and in- 
clinations (see Table [T]i. For our purposes below, however, 
it is useful also to have the ratio R\/a, where Ri is the ra- 
dius of the primary, and a the semi-major axis. This ratio 
is well constrained by the eclipse, depending on the eclipse 
duration 4 (at half maximum depth) and inclination / via 
sin^(7rfe/P)sin^/= (/;i/fl)2-cos2; (Russell 1912). We mea- 
sured tg/P graphically and list the values, as well as the in- 
ferred ratios Ri/ain Table[T] Note that the derivation assumes 
that the stars are spherical; if they are rotating rapidly, th en th e 
constraint is mostly on the equatorial radius (see Section lZSb . 

The measurement of R\/a, combined with Kepler's 
law, constrains the mean density of the primary, = 
(/?i/a)-^(37r/Gp2[l +q]). For both KOI 74 and KOI 8 1, the 
mass ratios q = M2 /Mi are small, and, as discussed by IRIOI 
the densities inferred for q ^ I, of ^^0.45 and 0.17gcm"^, 
respectively, are consistent with those expected for main- 
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Figure 1. Kepler light curves for KOI 74 (left) and KOI 81 (right). The top panels show the raw photometric data, and the bottom panels the detrended and 
folded light curves. The larger dips correspond to eclipses of the hot companion, while the smaller ones correspond to transits over its parent star. The smooth 
curve represents a fit to the blue data points - those taken outside of transits and eclipses, and within 3(t. The fit includes two sine waves, at the orbital and 
half-orbital period, representing Doppler boosting due to orbital motion and ellipsoidal light variations, respectively. 



sequence stars with masses of 2.2 and 2.1 Mq, respectively, 
as inferred from their spectral types (see Table[Tll. 

The ratio of the transit and eclipse depths et and eg equals 
the ratio of the surface brightnesses. Thus, given the tem- 
peratures for the primaries inferred from their spectra, one 
can estimate the temperatures of the secondaries. Assuming 
blackbody emission, one has. 



(1) 



Using the numbers from Table [T] and assuming an effective 
wavelength (A) = 6000 A, we infer T2 = 13000 and 17000K 
for KOI 74 and KOI 8 1 , respectively, with uncertainties some- 
where between 5% and 10%. These temperatures are higher 
than inferred by RIO, who assumed that the eclipses and tran- 
sits measured bolometric surface brightness. (We note that, 
of course, one could derive more precise values by folding 
blackbodies or model atmospheres through the Kepler pass- 
band, but this would not necessarily be more accurate as long 
as one does not account properly for, e.g., gravity darkening 
on what may well be a rapidly rotating primary. In any case, 
for our purposes the above estimates suffice.) 

2.2. Ellipsoidal Variations and Doppler boosting in KOI 74 

For KOI 74, in addition to the long-term trend, eclipses, and 
transits, the raw fluxes show sinusoidal modulations at the or- 
bital and half the orbital period (Figure [U. We fitted the data 
outside the eclipses and transits with a 6th-degree polynomial 
and sine waves at the orbital frequency and its harmonic (and 
iterated once, rejecting > 3a outliers). In the bottom panel of 
the figure, where fluxes normalized by the polynomial com- 
ponent of the fit are shown as a function of phase, the modula- 
tions are very obvious. They have much larger amplitude than 
is seen in Fig. 1 of RlOvl, likely because these authors origi- 
nally had removed most of the signal in their pre-display me- 
dian filtering; their later version shows a larger signal, more 
similar to what we find. (Related, we note that the two archive 
light curves include a 'corrected' flux. In the second, longer 
light curve, the signal is clearly present, but in the first, shorter 



one, it is not, apparently having been removed by the 'correc- 
tion' applied in the Kepler pipeline.) 

From our fit, we infer fractional amplitudes A\ = (1.082 ± 
0.013) X lO-"* and Ai = (1.426 ±0.018) x lO'^ for the fun- 
damental and the harmonic, respectively. As can be seen in 
Figure [T] these two fit the data quite well; although the fit is 
not formally acceptable (Xred =1-5 for 1843 degrees of free- 
dom), the residuals appear quite white. A Fourier transform 
of the residuals confirms this impression; there is no power 
above fractional amplitude 10"^, apart from one faint signal 
with ^ ~ 1.6 X 10"^ and / ~ 1.66d"'. In Figure[Tl one sees 
that the maximum of the fundamental coincides with one of 
the maxima of the harmonic; our fit gives a phase difference of 
-0.009 ± 0.003 cycles (of the orbital period; the maximum of 
the fundamental occurring slightly earlier). The fundamental 
has essentially the same offset (-0.008 ± 0.002 cycles) from 
the descending node of the massive star, while the harmonic 
has maxima coincident with the nodes to within the measure- 
ment errors (0.0005 ± 0.0008 cycles). 

The above precisions are remarkable, and hence we 
checked whether the error estimates were reliable. We found 
that the amplitudes do depend on how we model the long- 
term trend; if we reduce the polynomial to 4th degree, we find 
A12 = (1.034, 1.430) X 10""*; if we fit the pipeline 'corrected' 
fluxes in the longer light curve with a constant and the two 
sine waves, we obtain Ai,2 = (1-076, 1.515) x 10""*. This sug- 
gests that the true uncertainty in the amplitudes is about 5%, 
although we note that these latter fits are not as good (we find 
= 1.7, and the phases do not agree as well with expecta- 
tions). Obviously, the uncertainties will reduce substantially 
as further data are added. 

Both signals have natural if, for the fundamental, unusual 
interpretations. We believe that the fundamental is due to or- 
bital Doppler boosting, an effect which, as far as we are aware, 
was first detected and discussed by Maxted et al. (200 01); it is 
disc ussed in the c ontext of Kepler by Loeb & Gaudii (l2003h 
and IZucker et al.l (l2007h . Here, the precision of the Kepler 
data allows us to go beyond detection, and use the Doppler 
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boosting to measure a star's radial-velocity amplitude. 

Briefly, Doppler boosting occurs because as the primary or- 
bits, its spectrum is Doppler shifted, its photon emission rate 
is modulat ed, and its em i ssion is slightly beamed forward. 
Following Loeb & Gaudil (120031) . we use the relativistic in- 
variant quantity Ii,/v^ (or where I^, is the specific in- 
tensity, and V the frequency. Then, for a radial velocity v, , the 
photon rate observed by a telescope with given effective 
area Ax is given by, 

where F\ = Ix dft is the observed flux (fi is the solid angle 
subtended by the star; note that the integral is correct only to 
first order in v^/c). For a narrow band, the fractional variation 
would then be 

— — =/dB-- (^-5--^^j-~ ^hc/XkT_l 7' 

(3) 

where we introduced a Doppler boost pre-factor /db and 
where the approximate equality is for blackbody emission. 
For KOI 74 (r = 9400K), observed with Kepler ((A) ~ 
6000 A), one thus expects a signal with fractional amplitude 
of ^2.8(/r/c), where K is the radial-velocity amplitude. For a 
more precise estima te, we folded a 9500 K model atmosphere 
jMunari et al.ll2005h through the Kepler response. This yields 
/db = 2.25; the difference with the blackbody results from the 
fact that the model continuum is closer to Rayleigh Jeans in- 
side the Kepler bandpass, mimicking the emission from a hot- 
ter blackbody. In this temperature range, the boost pre-factor 
scales approximately linearly with temperature; for T = 9000 
and 10000 K, we find /db = 2.37 and 2.15, respectively. In- 
cluding reddening of E(B-V) = 0.15 (from the Kepler Input 
Catalog) has only a small effect, yielding /db = 2.21. Given 
that adding the extinction is equivalent to a relatively large 
change in effective area, the small change in /db found also 
means that uncertainties in the response are not important. 

Given the above, using fuB = 2.21, our observed light am- 
plitude corresponds to a radial-velocity amplitude of 14.7 ± 
1.0 kms"', where we assumed 5% uncertainties in both 
A\ and in /db- The corresponding mass function is AA = 
0.0017 ± O.OOO4M0. If we adopt a mass for the primary star 
of 2.22M0 (IRIO; see SectionlZTTi. thi s yiel ds M2 ~ 0.22Mo. 

Turning now to the harmonic, as in IRIOI we interpret it as 
changes in the observable surface area due to tidal distortion. 
The expected amplitude is given by, 

^2 = /EVTf - sin^-, (4) 
Ml \ a J 

where /ev is a pre-factor of order unity that depends on 
the limb darkening and gravity darkening. For our system, 
the eclip se and transit light curves yield / and Ri/a (Sec- 
tion l2.1b . Combining these with the amplitude yields a mass 
ratio M2/M1 = 0.092//ev- If we again adopt a pr imary mass 
of 2.22Mq, and take /ev — 1.5 (see Section |23] |. we derive 
M2 — O.14M0, somewhat below the value obtained above 
using the radial velocity amplitude, but roughly consistent, 
given the probably larger uncertainties associated with deter- 
mining the mass ratio via ellipsoidal light variations (see Sec- 
tion |22]l. 



2.3. Other Signals 

Apart from the signals discussed above, one might also ex- 
pect to see contributions of other factors, in particular the 
Doppler boosting and ellipsoidal variations of the secondary, 
irradiation effects on both primary and secondary, and possi- 
bly influences from an eccentric orbit. Following the proce- 
dure above, we find that the fractional amplitudes for the fun- 
damental and harmonic of the secondary, relative to the sec- 
ondary's flux, are ^ 10"^ and ~2 x 10"^, respectively (where 
we used K2 ~ 150 kms"', /db = 1.9, R2/a ~ 0.0026, and 
/ev = !)• Since the secondary's flux is only 0.2% of the to- 
tal flux, however, the contribution to the observed signal is 
negligible. 

For irradiation of the primary by the secondary, one expects 
a change in luminosity of ALi ~ TrR^iLa/ATra^), correspond- 
ing to a fractional amplitude jALi/L = ^(L2 / L){R\ /2a)^ ~ 
3.0 X 10"* (where L ^ L\ is the total luminosity). Within 
the Kepler passband, this should lead to a modulation of 
^ 1 .7 X 10~*. Similarly, for the secondary, one expects a frac- 
tional amplitude of i(7?2/2a)^ ~ 0.8 x 10"'' in bolometric flux, 

and a signal of '--^0.4 x 10~* in the Kepler passband. The two 
will have opposite phase; thus, the net expected signal has an 
amplitude of ^ 1.3 x 10"*, and is phased such that maximum 
occurs at the time of transit. This signal may be responsible 
for the fact that the fundamental is slightly out of phase with 
what is expected from Doppler boosting; the observed out- 
of-phase amplitude is ^1 sinA(/)i = (9±3) x 10"^, consistent 
with what we derived above. 

Finally, we consider possible signals from a non-circular 
orbit. For small eccentricity, both the Doppler boost curve 
and the ellipsoidal variations would no longer be pure si- 
nusoids, but have some contribution from higher harmonics. 
For the easier case of the Doppler boost curve, we can con- 
strain the eccentricity from the absence of a signal at the first 
harmonic that is out of phase with that expected from ellip- 
soidal variations (to a Hmit of A02 = 0.0005 ± 0.0008 cy- 
cles; Section |Z2] ). We find that, at 95% confidence, esinw = 
(^2/-4i)sin47rA</)2 < 0.03 (where lo is the longitude of peri- 
astron, and this expression is valid in the limit of small e). One 
sees that one cannot constrain the case where periastron oc- 
curs at one of the nodes, since in that case the harmonic would 
be hidden by the ellipsoidal variations. A complementary 
constraint, however, comes from the fact that the mid-transit 
and mid-eclipse times are consistent with a circular orbit. The 
measured offset from half a cycle is A0 = 0.0003 ± 0.0008 
cycles, which, at 95% confidence, implies ecosw = nAcf) < 
0.006 (with the expression again being valid in the limit of 
small e). Thus, from both constraints combined, we conclude 
that the orbit is circular to within e < 0.03. 

2.4. Application to KOI 81 

For KOI 81, a detailed analysis is not possible, since varia- 
tions on its much longer orbital period are strongly covariant 
with the trends in the raw data. Furthermore, the star is clearly 
pulsating; a Fourier transform shows that there is signal at nu- 
merous frequencies. In order to look for orbital modulation, 
we fit the raw fluxes with orbital and half-orbital modulations, 
a 3rd degree polynomial, as well as sine waves at the five 
largest pulsational signals, at 0.362, 0.723, 0.962, 1.32, and 
2.08 d"^ (some of these are harmonically related, as is often 
the case for multi-periodic pulsators). From the fit, orbital 
modulation does appear to be present, as is also clear in the 
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Table 1 

System parameters 



KOI 74 KOI 81 

Property Primary Secondary Primary Secondary 

Porb(d) 5.18875±0.00008 23.8776 ± 0.0020 

((°) 88.8 ±0.5 88.2 ±0.3 

Eclipse duration (cycle) 0.0362 ± 0.0004 0.058 ± 0.004 

Eclipse depth (10-5) 51^5 118 ±5 160 ±5 496 ±5 

R/a 0.116±0.002 0.0026 ± 0.0002 0.058±0.004 0.0023±0.0003 

Spectral type Al V ■■■ B9-A0V 

Radius (R0) l-90!^05 0.043±0.004 2.93±0.14 0.117±0.012 

7;ff(K) 9400 ±150 13000 ±1000 10000 ±150 17000 ±1300 

Luminosity (Lq) 25.6 ±2.4 0.05 ±0.02 77.3 ±9.6 0.9 ±0.4 

Velocity amplitude ( km s"' ) . . . 14.7 ± 1 .0 ■ ■ ■ ~ 7 

Mass(M0) 2.22!^-|5 0.22±0.03 "^-I^^W ~0-3 

Model mass (Mq) 0.20 ±0.03 ■■■ 0.25 ±0.03 

Note. — The periods and inclinations are taken from IRowe et al.l )2010i0 . as are the spectral types of the primaries 
and the quantities inferred from those. Systematic errors are possible if the star has accreted significant matter, or is 
rotating rapidly. The eclipse durations and depths were measured graphically, and have "chi-by-eye" uncertainties. For 
the temperatures of the secondaries, we assumed an 8% systematic uncertainty in our conversion of surface brightness 
to temperature. The velocity amplitude is measured from the Doppler boost signal, and its uncertainty includes possible 
systematic effects related to the detrending of the data and the precise spectrum of the primary. The model mass is based 
on the theoretical relation between orbital period and WD mass (see Equation js)). 

folded, normalized light curve, where we divided the fluxes Fitting the raw data of both sources to our model, we con- 
by the fitted trend and pulsational signals (see Figure[TJ- Our firmed that for KOI 81, the results depend sensitively on how 
fit yields amplitudes Axc^Sy. and ^2 ^ 4 x 10"^ which one fits the long-term trend. We thus decided to focus on 
corresponds to ~ Vkms"' and q ~ 0.2//ev. Because of KOI 74. We searched over a grid in q, Ku and Ti„, and fit- 
the long orbital period, the results are very sensitive to, e.g., ted for the best 4th-degree polynomial at each position. The 
the degree of the polynomial used. Comparing with a fit us- best-fit synthetic light curve is nearly indistinguishable from 
ing a 6th degree polynomial, as well as one using the 'cor- the two-sine fit, and has the parameters q = 0.0689 ±0.0009, 
rected' flux from the Kepler ai'chive light curve (the latter for K = 14.2 ± 0.2kms"', and ?]„ = 760 ± 40 K. We did not ex- 
the second, longer light curve only), the amplitude of the fun- elude outliers, and therefore our fit is poorer than the ana- 
damental changes by up to a factor of 2, while that of the lytical one (x^d - 2-6); we scaled the uncertainties such that 
harmonic changes by -20%. Furthermore, we find relatively ^2^^ = 1 (but did not attempt to include uncertainties in effec- 
large phase offsets, of ^-^0.03 cycles. Given this, we will not tive temperature, etc.). 

try to infer quantities, but simply note that for a 2.7M0 pri- These results reproduce our analytical estimates. The off- 

mary with a O.SMq companion, the predicted amplitudes are set in Ki would be even smaller if we had included reddening 

~8 X 10"^ and 3 x 10"^, respectively, consistent with the ob- in our numerical model (as this affects /db; see Section |2!2]|. 

servations. The mass ratio matches our estimates for /ev = 1.34, quite 

close to the value of /ev = 1-63 inferred from the tables of 

2.5. Verification Using a Light Curve Synthesis Code rflflg ). The secondary luminosity inferred from the ir- 

We tried to verify our semi-analytical estimates above us- radiation is 0.08 ±O.O2L0 (for orbital separation a = 16ARq 

ing a light-curve synthesis code, which is similar to that of and A = 0, i.e., assuming complete absorption and reradia- 

lOrosz & Hauschildfl (12000 1). and has b een used previously t o tion), consistent with our estimate from the surface brightness 

model irradiated pulsar companions (iStappers et alj 119991) . ratio. Below, we give our best estimates of the masses based 

The code accounts for tidal distortion and stellar rotation on these results, but first we discuss possible uncertainties due 

in the Roche approximation, includes irradiation and grav- to rapid rotation. 

ity darkening, and calculates the flux using NextGen model tt ^ ■ ^- j ^ n ^ ^■ 

; u JTJ ....^ .11 — -.or, irTnnn^ • . . j 2.6. Uncertainties due to Rotation 

atmospheres (IHauschildt. AUard. & Baronlll999.) , mtegrated 

over the Kepler passband. It cannot yet deal with eclipses. In many ways, the light curves of KOI 74 and KOI 81 are 

and hence we only attempt to reproduce the orbital modula- easy to model, since all effects are smaU and can thus ade- 

tions. quately be treated as perturbations. The one parameter that 

As input parameters, we chose to use the set Poib, kh U Tu we have not yet constrained, however, is the rotation of the 

Ri/a, Ki , f i.rot, q, and Tin, where the ones not mentioned be- primai-y. Since rapid rotation is expected in the context of the 

fore are to, the time of conjunction, and Ti„, the effective tem- evolutionary scenario described below, we discuss the obser- 

perature corresponding to the secondary flux absorbed by the vational consequences in some detail. 

primai-y (i.e., aT^ = (1 -A)L2/47ra^, with the albedo A ex- Rotation of the primai-y has a number of effects. First, the 
pected to be close to zero for a radiative star). Our main goal constraints on inclination and radius of the primary from the 
is to constrain q, K, and T,„; we assume the other parameters eclipse and transit change: the true incHnation /' will be closer 
are as inferred from the eclipses. For our fits, we also assume to 90° than the inclination / inferred assuming a spherical star, 
that the star is corotating with the orbit (Pi = Po,h); we will the radius will be smaller, and the radius that is constrained is 
return to this in SectionES] We add to the model the tiny flux the equatorial radius R^^ (assuming aligned rotation). Fur- 
contribution from the secondary (taken to be constant at the thermore, the estimate of the mean density increases, to p\ ~ 
level indicated by the eclipses). (^eq/a)~^(37r/GP^[l +^])(7?eq/-Rpoie), where Rpoie is the po- 



Doppler boosting in Kepler light curves 



5 



lar radius; thus, one would infer a somewhat smaller primary 
mass. For KOI 74, which has an inclination very close to 90°, 
we conclude that it would have /?eq/fl ^ Ri/a. For KOI 81, 
7?eq/fl would be somewhat smaller than Ri/a inferred from 
the eclipses. 

The second effect of a rapidly rotating primary is that its 
pole will be hotter than the equator; hence, looking at the 
equator, one will underestimate its temperature, and thus its 
luminosity and mass. This also affects the estimates for 
the secondary, since the ratios of the radii {R2/R1), based 
on the transit depth, and that of the surface brightnesses 
(Equation ([TJ) will be different from their true values: the 
cooler equator will lead to a larger inferred R2/R1 and smaller 
inferred 72. For the companion radius R2, however, the 
change in R2/R1 is partly compensated by the fact that a 
rapidly rotating primary also would have a smaller area {R[ ~ 

(^eq^poie)''^^; scc abovc). We have used our model code to 
test an illustrative case of extreme rotation, with the primary 
rotation rate 20 times faster than the orbital one (correspond- 
ing to vsin/ ~ 370 km s"'), and a ratio of equatorial to polar 
radius Rp, n/Rrm]p. = 1-32 (similar to what is measured for Reg- 
ulus; see iMcAlister et ani2005h . Keeping the flux- weighted 
temperature at 9400 K, we find equatorial and polar temper- 
atures of ^7800 and 11 600 K, respectively. For this particu- 
lar example, we estimate that the radius of the hot compan- 
ion would be R2 — 0.045 /?0, i.e., increased by only 5% from 
that inferred for a slowly rotating primary, and its temperature 
^ 10500K. 

Rapid rotation also indirectly affects the ellipsoidal varia- 
tions. First, because of the different temperature distribution 
on a rapidly rotating primary, different regions are weighted 
differently, and as a result the amplitude of the ellipsoidal 
variations changes. From our model code, we find that for 
the above-mentioned case, the predicted amplitude of the el- 
lipsoidal variations increases by a factor ^ 2. The effect is 
strongly non-linear, however; reducing the rotation rate from 
20 to 15 times faster than orbital, the amplitude increases only 
by a factor of 1.3. 

Another indirect effect of a rapidly rotating primary may 
be that the inferred q does not necessarily correspond to 
the true mass ratio. An implicit assumption underlying the 
equilibrium tide is that the primary corotates with the orbit. 
For stars that do not rotate sy nchronously, the work of, e.g., 
iPfahl. Arras. & PaxtonI (l2008h . suggests that the equilibrium 
tide may be a rather poor approximation to reality. Thus, it 
may be that one cannot reliably infer mass ratios from ellip- 
soidal variations for systems that are not tidally locked. 

We note that the light curves may contain a clue to the ro- 
tation rates. In particular, for KOI 74, we found evidence for 
a weak modulation at ^0.6 d, which is about 9 times faster 
than the orbital period, and corresponds to 30% of break-up 
for a 2.2Mq and \.9Rq star. For KOI 81, it may be possi- 
ble to infer the rotation rate from the pulsations by looking 
for rotational splitting (though for high rotation rates this is 
non-trivial). 

Of course, it would be simpler and more reliable to measure 
the rotational broadening spectroscopically, and the spectra 
mentioned by lRlOl may already hold the answer, but we do not 
have access to those. Nevertheless, we note in closing a much 
more subtle consequence of rapid rotation, which is t he pho- 
tometric analog of the Rossiter-McLaughlin effect (Rossite? 
[1924; McLaughlin 1924). Assuming the rotation is aligned 
with the orbit, at the start of the transit more blueshifted light 
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Figure 2. Representative evolutionary track in tlie Hertzsprung-Russell (HR) 
diagram, illustrating the evolution of the progenitors of the hot companions 
of KOI 74 and KOI 8 1 . The solid black curve shows the track of a star with 
an initial mass of 1.4M0 that starts to transfer mass to a companion star just 
after leaving the main sequence and ultimately becomes a 0.21 M0 helium 
WD. Because of residual hydrogen in the envelope at the end of the main 
mass-transfer phase, the degenerate remnant experiences three hydrogen shell 
flashes, producing three large loops in the HR diagram (of progressively 
larger area), before it settles on a classical He WD cooling track. Selected 
ages since the end of the main mass-transfer phase (labelled as "At = Oyr") 
are given next to filled circles. The dashed Unes give lines of constant radius 
as indicated. The positions of KOI 74b and KOI 81b are also shown with 
conservative error bars. They are consistent with the expected location for 
~ 0.2-0.25 Mq He WDs either on a classical He WD cooling track or during 
one of the relatively long-lived phases of a H-flash loop. 

is being blocked and at the end more redshifted light. For 
a rotational velocity of 300km s"', Doppler boosting induces 
a signal in the blocked light with a fractional amplitude of 
r^2x 10"^. This is diluted by the unblocked light, i.e., by 
a factor equal to the transit depth. Thus, for KOI 74 and 
KOI 81, one expects net signals of 1 x 10"^ and 3 x 10"^, 
respectively, which may be detectable by averaging about 100 
transits. 

2.7. Adopted Masses 

In principle, given our inferred values of Ki and q, com- 
bined with the inclination, one can derive both masses with- 
out further assumptions. In practice, however, these masses 
are very uncertain, since the uncertainties in Ki and q enter 
to high powers (for small q, Mi oc K^q~^ and M2 oc K^q~^). 
Given that these uncertainties are of order 5% for Ki and 
likely larger for q (see above), the la fractional errors are 
>20%, and the distributions are highly non-Gaussian. 

Instead, we proceeded by assuming the primary has a mass 
of 2.22^qJ^Mq, as inferred fro m it s spectral type and mean 
density jRlOl ; see also Section ITTl i. Then, using Ki and /, 
we calculate a mass of O.22±O.O3M0 for the secondary. If 
instead we use the mass ratio q, inferred from the ellipsoidal 
variations, we find M2 = O.15M0, with an uncertainty that is 
difficult to estimate because we cannot be sure the star is in 
corotation (see above). We will proceed by using the value 
inferred from Ki (see Table (TJ. For KOI 81, we do not have 
a good constraint on Ki (or q), but the signals that are present 
are consistent with a similar secondary mass, of ^0.3Mq. 

3. EVOLUTIONARY HISTORY 

An obvious explanation for these hot, compact companions 
is that they are WD stars. RIO have discussed this possibil- 
ity for KOI 81, but their initial very low mass estimate for 
the companion in KOI 74 (<O.O32M0; iRlOvl.) was inconsis- 
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tent with a WD explanation. Since our reanalysis of the data 
gives the mass of the companion in KOI 74 as ~ 0.2Mq, a 
WD co mpanion is no longer excluded. Cooling tracks from 
iPanei e t al. (2007) for helium-core WDs of appropriate mass 
are broadly consistent with the companion luminosities and 
temperatures in Table [T] (somewhat more easily for KOI 74 
than for KOI 81). Figure |2] compares a cooling track from 
our own calculation^ for a ~O.2M0 helium-core WD to the 
WDs in KOI 74 and 81. 

Binary systems resembling these may well be common; the 
well-known star Regulus (a Leonis) has recently been found 
to be a spectroscopic binary with componen t masses "^0.3 an d 
3.4Mq, and an orbital period of 40.11 d dGies et al.ll2"008h . 
The lower-mass component likely is a WD, making Regulus 
and its companion remarkably similar to the pair of systems 
considered here. The evolution of Regulus has been consid- 
ered in detail by lRappaport. Podsiadlowski. & HorevI {200^. 

To estimate the occurrence more quantitatively, we note that 
Kepler observ ed only ~400 stars with temperatures above 
about 9000 K (i Batalha et akl 120101) . Since the eclipse prob- 
abilities for systems like KOI 74 and KOI 81 are ~ 12% and 
6%, respectively, the detection of two such WD companions 
suggests that a surprisingly large fraction of such stars must 
have WD companions. With only two detections, it is difficult 
to estimate the orbital period distribution. Taking an average 
eclipse probability of ~ 9%, then for a fraction of A and 
B stars with suitably close WD companions, the mean num- 
ber of expected eclipses among the Kepler sample would be 
^ AQQFV. The Poisson probability for finding two or more 
eclipses out of 400 systems is ^50% for T = 0.05 and ^-^25% 
for F = 0.03. Therefore, we conclude that 1 out of every 20 
to 30 A and B stars is in a close, few tens of days, binary with 
a WD companion. 

The above abundance is very high, and, if confirmed, has 
significant implications for our un derstanding of stellar and 
binary evolution (see also iDi Stefa no 2010). For instance, 
the descendants of these systems should also be relatively 
numerous — and likely interesting — stars. When the currently 
more massive star leaves the main sequence and expands 
to overfill its Roche lobe, the resulting mass transfer will 
be dynamically unstable, and likely lead to a merger (even 
more likely than for R egulus, which has a longer period; see 
iRappaport et al.ll2009h . In this merger, the addition of a he- 
lium WD to the helium core of the subgiant star could lead 
directly to core helium ignition, bypassing the usual grad- 
ual buildup of the helium core via hydrogen shell burning. 
The merger product would also be rapidly rotating, and hence 
the natural outcome would appear to be a rapidly-rotating red 
clump or horizo ntal branch sta r (such as observed by, e.g., 
iBehr et al.ll2000l and lBehdl2003l) . 

Of course, while the above is interesting, it is possible that 
we are dealing with a statistical fluke. Indeed, if WD compan- 
ions are very common among A stars, it seems odd that none 
has been reported so far for any of the much more numerous 
cooler stars in the Kepler sample. On the other hand, this may 
indicate that binary systems with initially lower-mass stars, 
which would evolve to G and F stars with low-mass WD com- 
panions, have a difficult time forming and/or surviving. 

^ In the binary calculation shown in Figure |2] (from the library of 
IPodsiadlowski. R appaport. & Pfahl 2002), the companion is a neutron star. 
The tracks for the KOI 74 and KOI 81 systems are expected to be very 
similar (see also Rappaport et al. 200^. For more details on the hy- 
drogen shell flashes in these tracks, see IPodsiadlowski et ali 120021) and 
[Nelson, Dubeau, & MacCannell (200|). 
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Figure 3. Example evolutionary histories in the mass-peiiod plane. The red 
curves are for the mass-gaining stars, the blue curves for the mass-losing 
stars. The orbital periods of KOI 74 and KOI 81 are marked with green hnes. 
The M^j-Porb relation, where the mass transfer terminates, is marked by a 
black curve. Many different such potential histories could be constructed. 

Turning now to the formation of KOI 8 1 and KOI 74, for 
both the initial masses of the WD progenitors must have 
been greater than the initial masses of the current A stars, 
since the progenitors of the WDs have evolved more rapidly. 
Some form of binary interaction is needed in order to remove 
mass and produce the WDs. Either stable mass-transfer from 
the more massive to the less massive star has occurred; in 
this case, the current mass of the A star may be su bstan- 
tially higher than its initial mass (see, e.g., Rappapor t et alJ 
l2009h . Alternatively, the system may hav e passed through a 
common-envelope phase dPacz vnskill 1 976h . 

3.1. Stable Mass Transfer 

If the mass loss from the progenitor of the WDs were sim- 
ply through stable Roche-lobe overflow, then we could as- 
sume that the current orbital periods are very close to the 
orbital periods at the end of mass transfer, i.e., the point at 
which the mass of the envelope about each current WD be- 
came too small to maintain a giant structure and collapsedQ 
Hence we can apply the relationship between core mass and 
radius to predict the mass of each WD based on the cur r ent or- 
bital periods. Using this relationship. IRappaport et al] ( Il995h 
approximated the orbital period Porb as, 

Porb« 1.3x10^^^^^^^ days, (5) 

where the WD mass M„d is in units of solar masses. At a 
given Poib, the spread in M^d is expected to be at most ±15%. 

Applying Equation (|5]l we find M„d = 0.20 ± 0.03 M© for 
KOI 74, which matches the value independently derived in 
Section |2] For KOI 81, this becomes M„d = 0.25 ± 0.03 M©, 
again consistent with a low-mass WD0 

Stable mass transfer seems to provide a simple way to pro- 
duce both these systems. The match between the WD masses 
from Equation (|5]l) and from the light-curve data is encour- 
aging. Furthermore, the full binary evolution calculations of 

^ A-type stars typically experience no magnetic braking, which might oth- 
erwise have shortened the orbital period since the end of mass-transfer. 

Tauris & Savonije ( 1999) provide an alternative expression, from which 
we obtain similar values: M„d = 0.24 ± 0.02Mq for KOI 74 and M„d = 
0.29 ± O.O2M0 for KOI 81. 
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iPodsiadlowski. Rappaport. & Pfahll (l2002h produce systems 
with very similar orbital periods and WD masses (see their 
Fig. 13). In Figure |3] we show examples of how the two sys- 
tems could have reached their current state, us ing simple ana- 
lytic formulae for orbital evolution (Eq. (3) of Rappapo rt et alJ 
12009; see also Podsiadlowski, Joss, & Hsu 1992). The curves 
in Figure [3] take initial conditions {Mi init, Mji^ix, Porbinit} of 
{I.6M0, 1.2M0, 19.3h} for KOI 74 and {l.SM©, 1.3M0, 
31 h} for KOI 81. In both cases we assume that any mass 
lost from the system carries away the specific orbital angular 
momentum of the binary. For KOI 74, this example assumes 
that the mass transfer is 73% conservative, and for KOI 81 it 
is 91% conservative. There are certainly other possible paths 
to the current systems, in particular since the above ignores 
the possible role of magnetic braking (in evolutionary stages 
where either star has a convective envelope). 

The currently more massive stars are almost certain to have 
accreted significant amounts of matte r and are therefor e likely 
to be very rapidly rotating (see, e.g., iRappaport et al.l (120091) 
and references therein). Their spin angular momentum could, 
in principle, have been lost by magnetic braking. However, 
typical A-type stars do not experience strong magnetic brak- 
ing (see, e.g., Kawaler 1988). Alternatively, for KOI 74, tides 
might have been strong enough to slow the rotation of the star 
(and slightly widen the orbit). 

3.2. Alternative: Common-envelope Evolution? 

Common envelope evolution seems to provide a much less 
satisfactory explanation for these systems than stable mass 
transfer If the envelopes of the initially more massive stars 
were ejected during common envelope evolution, then us- 
ing the current Poib in Equation (|5]l would only yield a min- 
imum WD mass. So, if the observational data cannot allow 
more massive WDs than ~0.25Mq, it seems unlikely that ei- 
ther of these systems have experienced a significant common- 
envelope phase. 

In addition, any scenario in which the currently more mas- 
sive stars have not accreted a significant amount of mass 
seems difficult. The > 2Mq current primary stars in these 
systems would produce WDs more massive than >O.28M0 
if their envelopes were removed at the end of the main se- 
quence. If those initially less massive stars had not gained 
matter then, at best, significant fine-tuning would be required 
for the initially more massive star to have both evolved off the 
main sequence and produced the low-mass WDs observed. 

3.3. Other Expected Compact Hot Companions: Hot 
Subdwarfs 

Another class of hot compact companion stars that should 
result in light curves similar to those in Figure [T] are hot 
subdwarf (sdB, sdO) stars, sometimes referr ed to as extreme 
horiz ontal branch stars. The calculations of Han et al.l (120021 
120031) predict that sdB stars should exist in binaries of these 
orbital periods about A- and B-type stars (see lHan et al.ll2003l 
Fig. 15). However, at the temperatures inferred for KOI 74 
and 8 1 , sdB stars would be expected to be considerably more 
luminous. Typical hot subdwarfs properties are anywhere be- 
tween -20000 and 40000 K and -3 and lOOL©. In ad- 
dition, the inferred companion masses in KOI 74 is some- 
what too low to have ignited helium, even non-degenerately; 
masses greater than '-^O.3M0 are required. Such hot subd- 
warfs should be found in current and future planetary transit 
searches. 



4. SUMMARY 

The wonderful photometric precision of Kepler has al- 
lowed us to measure the radial velocity amplitude of KOI 74. 
Combined with the primary mass inferred from its spectral 
type and mean density, we estimated a companion mass of 
O.22M0, and argued it was a WD. We showed that its prop- 
erties are in very good agreement with theoretical expecta- 
tions based on a stable-mass transfer phase inside a binary 
system. A similar evolutionary history is also likely to have 
formed the current KOI 81 system. The many binaries that 
will be discovered by photometric surveys should be interest- 
ing and useful, both individually and collectively (see, e.g., 
IWillems et al.1120 06^ for a quantitative study), in helping us to 
increase our understanding of binary evolution. 

Future observations by Kepler should allow measurements 
of the photometric Doppler effect in more systems, perhaps 
including systems such as 'ultra-cool white dwarfs,' where the 
lack of spectral lines has precluded searches for radial veloc- 
ity variations. Given long enough baselines, the photometric 
analog of the Rossiter-McLaughhn effect also seems likely to 
be observed. 

Coffee-time astro-ph discussions at the Kavli Institute for 
Astronomy & Astrophysics triggered this work; in partic- 
ular, we thank Yanqin Wu, Andrew Shannon & Matthias 
Gritschneder. We thank the referee, Scott Gaudi, for his 
careful reading and useful comments. S.A.R. and Ph. P. thank 
J. Rowe for sharing his work in advance of submission. S.J. 
thanks Bill Paxton for the plotting package Tioga. 

Facilities: Kepler 

Note added in proof. After acceptance of this article, high- 
resolution spectra were taken for us by E. Kirby using HIRES 
on the Keck I telescope. From these, we measure projected 
rotational velocities of 150 and 225 km s' for KOI 74 and 
KOI 81, respectively, and infer rotational periods of about 
0.6 days for both sources. For KOI 74, the period is consis- 
tent with that inferred from the weak non- orbital modulation 
shown in the power spectrum (Section l272] i. 
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